Vortices, antivortices and superfluid shells separating Mott-insulating regions 
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Atomic or molecular bosons in harmonically confined optical lattices exhibit a wedding cake 
structure consisting of insulating (Mott) shells. It is shown that superfluid regions emerge between 
Mott shells as a result of fluctuations due to flnite hopping. It is found that the order parameter 
equation in the superfluid regions is not of the Gross-Pitaeviskii type except near the insulator 
to superfluid boundaries. The excitation spectra in the Mott and superfluid regions are obtained, 
and it is shown that the superfluid shells posses low energy sound modes with spatially dependent 
sound velocity described by a local index of refraction directly related to the local superfluid density. 
Lastly, the Berezinskii-Kosterlitz-Thouless transition and vortex-antivortex pairs are discussed in 
thin (wide) superfluid shells (rings) limited by three (two) dimensional Mott regions. 

PACS numbers: 03.75.Hh, 03.75.Kk, 03.75 Lm 
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Introduction: The recent experimental discovery of 
Bose-Mott insulating phases in optical lattices has gener- 
ated an explosion of research in the ultra-cold atom com- 
munity (see [H for a recent review), and has helped to 
merge two major branches of physics: atomic-molecular- 
optical and condensed matter physics. Most experiments 
thus far have relied on time of flight measurements to 
infer the existence of a superfluid to insulator transi- 
tion 0, JsL However, very recently, two experimental 
groups [j,l5| have used spatially selective microwave spec- 
troscopy to probe in situ the superfiuid-to- insulator tran- 
sition of ^^Rb in a three dimensional (3D) optical lattice 
with a harmonic envelope. In these experiments, the shell 
structure of the Bose-Mott insulating states was revealed 
for very deep lattices. Regions of filling fraction n — 1 
through n = b {n = 1 through n = 3) were mapped in 
the MIT [^] (Mainz [5(|) experiment. Their observation 
in three-dimensional optical lattices lead to the confirma- 
tion of the Mott-insulating shell structure consisting of 
"Mott plateaus" with abrupt transitions between any two 
sucessive shells as proposed in two dimensional (2D) Q 
optical lattices. 

One of the next frontier for ultra-cold bosons in opti- 
cal lattices is the search for superfluid regions separating 
Mott-insulating shells. These additional superfluid re- 
gions can be determined through the use of microwave 
spectroscopy much in the same way as the Mott shell 
structure was determined recently [j, 0|. Thus, in antici- 
pation of the next experimental breakthrough, we study 
2D and 3D optical lattices of atomic or molecular bosons 
in harmonically confining potentials, and show that be- 
tween the Mott regions of filling fraction n and n -f 1, 
superfluid shells emerge as a result of fluctuations due to 
finite hopping. This finite hopping breaks the local en- 
ergy degeneracy of neighboring Mott-shells, determines 
the size of the superfluid regions as shown in Fig [U and 
is responsible for the low energy (sound) and vortex ex- 



citations. In addition, we find that the order parameter 
equation is not in general of the Gross-Pitaeviskii type. 
Furthermore, in 3D optical lattices, when superfluid re- 
gions are thin (nearly 2D) spherical (or ellipsoidal) shells, 
we obtain bound vortex-antivortex excitations below the 
Berezinski-Kosterlitz-Thouless (BKT) transition temper- 
ature I?!, Isi] which is different for each superfluid shell. 
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FIG. 1: (color online) a) Shell structure of Mott and super- 
fluids regions in a 2D square optical lattice with harmonic 
envelope as a function of radius r/a for t 7^ 0. The super- 
fluid regions are shown in red (gray) whereas the Mott re- 
gions are shown in white. The black circles indicate the Mott 
boundaries Rc,n at t = 0. b) The local filling factor n(r) 
is shown in solid (dashed) lines for t ^ Q {t = Q). The red 
curve (solid gray) shows the local superfiuid order parameter 
|V'(r)p. The parameters are = 6 x 10"''?7, t = 1.25 x 10"^(7 
and ^ — 2.5U. 

Hamiltonian: To describe the physics of alternating 
insulating and superfluid shells for atomic or molecu- 
lar bosons we use the lattice Bose-Hubbard Hamiltonian 
with a harmonic potential described by 
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where j^i^ = 1^ — ^C*") is the local chemical potential, 
V{r) = Vlp{p/a)'^/2 + SX^X^jaf is the harmonically 
confining potential, a is the lattice spacing, and is 
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the creation operator for an atom at site r. Here, t is the 
hopping parameter and U is the interaction strength. 

We introduce a local superfluid order parameter ip^ = 
(cr) that produces an effective local Hamiltonian 

Hf - Ho,n (r) - < ^ (CrV^;+a + 4V^r+a ~ ^r>r+a) , (2) 



which is diagonal in the site index r with i7o,Ti(r) = 
^"hrifir ~ 1) ~ fJ-rf^r, whcrc fi^ — cjcr is the number 
operator. For t ^ 0, the Mott shell structure is re- 
vealed by fixing fir = n, to obtain the local energy 
Eo,n{^) — Un{n — l)/2 — /ir^i, when (n— 1)U < fir < nil. 
Since i?o,n+i(i") — £'o,n(r) — nil — fir, the change from a 
Mott shell with filling fraction n to n -I- 1 occurs at the 
degeneracy condition fir = nlJ , which for a spherically 
symmetric potential happens at i?c,„ = a^2{fi — nU)/n. 
The relation fir = nU determines the shape and size of 
the boundary between the n and n + 1 shells. However, as 
shown next, near this region of degeneracy fluctuations 
due to hopping introduce superfluid shells. 

Nearly Degenerate Perturbation Theory: We focus our 
attention now on the Mott regions with integer boson 
filling n and n + 1 and the superfluid shell that emerges 
between them. In the limit where U ^ t we can restrict 
our Hilbert space to number-basis states \n) and \n + 
1) at each site. Any contribution of other states to the 
local energy will be of the order of t^ /U . The hopping 
term in Eq. ([2]) affects the ground state energy of the 
system by removing the local degeneracy of -Eo,ra+i(r) 
and i?o,n(r) at fLr — nU. To illustrate this point, we 
make a continuum approximation and perform a Taylor 
expansion of ^'(r-l-a) to second order in a, provided that 
r ^ a. This approximation is not really necessary, but 
allows substantial analytical insight leading to 



jjeS^ ( Eo,r,{r) + A{r) 



-V?^H~TA(r) 

-V^7TTA*(r) £;o,„+i(r)+A(r) 
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where A(r) = i(A(r)V'*(r) -I- cc) and A(r) = t{z->p{r) + 
a^V^'(/'(r))- Here, z is the coordination number which de- 
pends on dimension d. Notice that t has two effects. 
First, it changes the local energies £'o,n(r) and £'o,ra+i(r) 
of the Mott shells n and n + 1 through A(r). Second, 
it mixes the two Mott regions through the off-diagonal 
term ^/n + lA(r) and its hermitian conjugate. Thus, the 
physics near the boundary between the n and n+1 Mott 
regions is described by an effective local two-level sys- 
tem with diagonal (A(r)) and off-diagonal {Vn + lA(r)) 
perturbations. 

The eigenvalues of Eq. ^ are given by, 

E^ir) = E^ir) ± ^ {E,{r)y + {n + 1) \Air)\\ 

where Es{r) = [£'o,„+i(r) + Ea^i^)] /2 + A(r) is pro- 
portional to the sum of the diagonal terms, and Ed{r) — 



[£'o,„+i(r) - -Bo,„(r)] /2 = {nil - fir)/2 is proportional 
to their difference. Notice that -E'-(r) is the lowest 
local energy, leading to the ground state energy E = 
^JdvE.ir). 

Order Parameter and Compressibility: The order pa- 
rameter equation (OPE) is determined by minimization 
E with respect to '0*(r) leading to 



A(r)- 



(n -h l)t(z + a^V^)A(r) 



= 0. 



(4) 



(n + l)|A(r)r 



Notice that the OPE is not of the Gross-Pitaeviskii (GP) 
type, since the superfluid regions emerge from local fluc- 
tuations between neighboring Mott shells. The zeroth 
order solution of this equation leads to the spatially de- 
pendent order parameter 
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Since |'0(r)| > 0, hence \nU — fir] < {n + l)zt, and 
the saddle point inner R„^- and outer Rn,+ radii for the 
superfluid shell between the n and n + 1 Mott regions is 
obtained by setting |V'(r)|^ = leading to 



Rn.± = Rc.n\ 1 ± 



2zt{n+l) 



n 
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for a spherically symmetric harmonic potential V{r) — 
r2(r/a)^/2, where i?c,n is defined above. This rela- 
tion shows explicitly that t splits the spatial degener- 
acy of the n and n + 1 Mott shells at r = i?c,n or 
{fLr = nil) by introducing a superfluid region of thick- 
ness Ai?„ = Rn,+ — Rn.-- In the case of non-spherical 
harmonic potential V{r) = Qp{p/a)^/2 + ^z{z/a)'^/2 the 
shell regions are ellipsoidal instead of spherical. Notice 
that ARn depends strongly on filling fraction n, chemical 
potential fi through Rc,n and the ratio zt/fl. 

In addition, the local filling fraction 



n{r) = - 



dfi 



= n+ - - 



nil — fir 
2zt{n + l) 



(7) 



in the same region interpolates between n+1 for r ^ Rn.- 
and n for r ^ Rn,+, while the chemical potential fi is 
fixed by the total number of particles N = J drn{r). In 
Fig. m n{r), \i^{r)\\ and i?„^± are shown for the Mott 
and superfluid regions. The local compressibility 



K(r) 



dn{r) _ 1 
dfi ~ 2zt{n+ 1) 



(8) 



of the superfluid shells is non-zero, in contrast to the in- 
compressible (k = 0) n and n+1 Mott shells for r < 
and r > Rn.+ , respectively. The local compressibility in- 
dicates the presence of large (small) excitation energies 
in the Mott (superfluid) regions, as discussed later. 
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Gross- Pitaeviskii behavior at the edges: We note in 
passing that only near the edges of the superfiuid region 
(where r « Rn,± a-nd "0(1") ~ 0) a direct expansion of 
the OPE (Eq. 2]) leads to the effective Gross-Pitaeviskii 
equation 



V2 + yeff(r)+geff|i/'(r)rU(r)=0. (9) 



2moff 



Here, rrics = l/ia^t is exactly the boson band mass, 
14ff(r) = \nU — ^r|/(?i + 1) — zt is the effective poten- 
tial, and Qcs — 2zt/{n + 1) is the effective interaction. 
Notice that Voff(r) < and vanishes at the boundaries 
Rn,± of the superfiuid region since \nU — /ir| = zt{n + 1) 
there. Furthermore, notice g^s = zt/{n + 1) is small in 
comparison to U , indicating that the superfiuid near the 
edges is weakly interacting, and more so as the Mott 
index n increases. When i — > (moff — > oo) then 
|'0(r)P ~ —Vcs/gcS leading to the correct limiting be- 
havior of Eq. m near r w Rn.±- Next, we discuss the 
excitation spectrum in the Mott and superfiuid regions. 

Excitation Spectrum: The excitation spectrum in the 
Mott and superfiuid shells can be obtained using the 
functional integration method 9] for the action {h = 
fcs = 1, /3 = l/T) 



leading to 
/ Pc^Pcexp - 



,t ^1 _ 



/3 



dr^ [cl,drC,,r + H] 



(10) 



the 
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partition function Z = 
In each Mott shell (away 
from the superfiuid region) we introduce a Hubbard- 
Stratanovich field ^' to take into account fluctuations 
due to the presence of finite hopping, and integrate out 
the bosons {c\ c) leading to an effective action 



{iuj, r) 



to quadratic order in and Vf, where k, k' are momen- 
tum labels, uj are Matsubara frequencies, and 



with -Ei(r) 



1 



nU 



ioj — Ei(r) 
Hr, E2{r) = {n 



E2{r) 



1)U — iir, and 
ek — 2tJ2'j=i'^os{kja). The poles of G'kk'(*^, r) are 
found upon the analytical continuation iuj = uj + i6, lead- 
ing to the local excitation energies 



-Mr 



+ - (2n - 1) - 1 ± _ (4n + 2)6kt/ + C/2 



where the +(— ) sign labels quasiparticle (quasihole) ex- 
citations. The energy to add a quasiparticle is Eqp — lu^, 
while the energy to add a quasihole is Eqh — —i^-- 

Figure ([2]) a shows the quasiparticle and quasihole ener- 
gies of the Mott phase as a function of position r. The en- 
ergy cost to create a quasiparticle (quasihole) is minimum 
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FIG. 2: (color online) a) Quasiparticle Eqp (solid line), quasi- 
hole Eqh (dashed line), and Mott gap Eg (red) energies for 
k = versus r/a b) sound velocity for the outermost super- 
fluid ring versus r/a. Same parameters as in Fig. [1] 



(maximum) at the trap center and increases (decreases) 
radially, while the Mott- Hubbard gap Eg = Eqp + Eqh is 
large and independent of r. 

The excitation spectrum in the superfiuid region can 
be calculated by a similar method. First, we intro- 
duce the Hubbard-Stratanovich field ip which now cor- 
responds to the order parameter in the superfiuid re- 
gion, second we use an amplitude-phase representation 
^{r, t) — |7/'(r, r)| exp [iip{v, r)] and we apply the nearly 
degenerate perturbation theory described earlier to inte- 
grate out the boson fields and c. Thus, we obtain the 
phase-only effective action 



'S'e 



1 

21? 



dvdr [k [dr^Y + Pij y^^j ^p] (11) 



to quadratic order in the phase variable for the superfiuid 
region between the n and n+1 Mott shells. Here, we as- 
sumed that |'0(r, t)| is r-independent at the saddle point. 
The coefficient k is the compressibility of the superfiuid 
described in Eq. O and 



p,j = 2ta^\ij{r)\^Sij = ps{r)5r 



(12) 



is the local superfiuid density tensor [lo| . which is di- 
agonal and vanishes at the Mott boundaries i?„.±- The 
resulting wave equation has the form 



0, 



(13) 



leading to a local sound velocity c(r) — ^/pjjjjfi, 
which in terms of the order parameter reads c(r) = 
2^/(rl+T)zta\^^jJ{'r)\. The local speed of sound has it max- 
imal value at |^/'(r)|m(ja; = (V" + l)/2, and the superfiuid 
region behaves as a medium of continuous index of refrac- 
tion x(r) = Cmax/c{r) = + l/(2|i/'(r)|). Notice that 
x(r) ^ oo at the Mott boundaries where |V'(r)| = 0, in- 
dicating that the sound waves of the superfiuid do not 
propagate into the Mott regions. A plot of the local 
sound velocity is shown in Fig. for the superfiuid re- 
gion between the n = 1 and n = Mott shells. From 
the phase-only effective action for the superfiuid region, 
we can also investigate the vortex excitations, which are 
discussed next. 
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Vortices: Next, we explore vortex solutions in two 
cases where spontaneous vortex-antivortex pairs can ap- 
pear as indicators of the Berezinski-Kosterlitz-Thouless 
(BKT) transition 0, HI. Case I corresponds to a 3D sys- 
tem, where the superfluid regions are very thin A_R„ <^ 
Rem leading to essentially a two-dimensional superfluid 
in curved space. Case II corresponds to a 2D system, 
where the superfluid regions are thick rings Ai?„ Rem 
leading to essentially a two-dimensional superfluid sub- 
ject to boundary conditions. 

In a flat space two-dimensional system stationary vor- 
tex solutions must satisfy §\Jip-d\ — 27rm, where 
TO = ±1,±2, ... is the vorticity (topological charge) and 
V(/9 is the superfluid velocity. The standard vortex so- 
lution in cylindrical coordinates is V(/3 = ni9/r, and 
the corresponding free energy per unit volume is T — 
2^ / drps (r) (Viy9)^ . This situation is analogous to a two- 
dimensional linear dielectric material where the displace- 
ment field is D = Vtpxz = e(r)E, with dielectric function 
e(r) = 1/ [27rp(r)]. Notice that the dieletric function di- 
verges at superfluid boundaries, since ps(r) — > there. 
In this language T is identical to the electrostatic energy 
per unit volume U^y — I ' general the so- 

lutions for several vortices (antivortices) can be obtained 
from VAV(/? — 27rz; TOi(5(r — r^), where is the location 
of the vortex (or antivortex) of vorticity rrii. 

In case I the superfluid state appears below Tbkt ~ 
Trps{r = Rc.n)/2, where ps = Ps/ol^ has dimensions of 
energy. In this limit ps{r = Ren) — {n + l)t/2 and the 
critical temperature of the superfluid shell between the n 
and n + 1 Mott regions depends on the index n. Notice, 
however that such estimate is reasonable only when the 
shells are sufficiently large. The solution for a vortex- 
antivortex (VA) pair in curved two-dimensional space is 



(p{r = Rc,n,0, 



arctan 



46i?c,„tan (^) sin((?!)) 
^ 62_4^2^tan2(^) 



where Rc.n = a\/'2{iJ. — nU)/il is the radius of the su- 
perfluid shell, and b is the VA size. A three-dimensional 
view of the velocity field V(p{9, 0) is shown in Fig. [S] 
When the superfluid shell has a small thickness Ai?„ then 
Tbkt ~ Trt{n + l)Ai?„/6a, while the vortex-antivortex 
pair has an approximate solution of the same form as 
above which interpolates between ip{r = Rn_^,9,(j)) and 
ip{r = Rn,+,9,(p). 

In case II the superfluid regions are rings bounded by 
Rn.- and Rn.+ , and one can use the Coulomb gas anal- 
ogy described above, conformal mapping techniques and 
proper boundary conditions to obtain vortex-antivortex 
solutions. The creation of vortex-antivortex pairs are en- 
ergetically quite costly when A_R„ ^ Rc.n, due to strong 
confinement effect of the boundaries, thus we do not ex- 
pect a BKT-type superfluid transition to occur until Ai?„ 
is substantially large (~ Rc,n)- Only in this limit, we 
expect a BKT transition with Tbkt ~ 7r(A5)/2, where 



{ps)/2 is the surface area average of ps{r) 



FIG. 3: Three-dimensional view of a vortex-antivortex pair 
in 2D superfluid shell separating two 3D Mott regions. 



Conclusions: We studied 2D and 3D optical lattices 
of atomic or molecular bosons in harmonically confin- 
ing potentials, and showed that between the Mott re- 
gions of filling fraction n and n + I, superfluid shells 
emerge as a result of fluctuations due to finite hop- 
ping. We found that the presence of finite hopping 
breaks the local energy degeneracy of neighboring Mott- 
shells, determines the size of the superfluid regions as 
shown in Fig[l] and is responsible for low energy (sound) 
and vortex excitations. In addition, we obtained an 
order parameter equation which is not in general of 
the Gross-Pitaeviskii type. Furthermore, we obtained 
bound vortex-antivortex solutions (as shown in Fig[3]be- 
low the Berezinski-Kosterlitz-Thouless (BKT) transition 
when superfluid regions are thin (nearly 2D) spherical 
(or ellipsoidal) shells. The emergence of these super- 
fluid regions should be detectable experimentally using 
microwave spectroscopy 0, Q . 
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